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Recent results of source function imaging from AGS through CERN SPS to RHIC 
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Recent femtoscopic measurements involving the use of an imaging technique and a newly developed moment 
analysis are presented and discussed. We show that this new paradigm allows robust investigation of reaction 
dynamics for which the sound speed c s ^ during an extended hadronization period. Source functions extracted 
for charged pions produced in Au+Au and Pb+Pb collisions show non-Gaussian tails for a broad selection of 
collision energies. The ratio of the RMS radii of these source functions in the out and side directions are found 
to be greater than 1, suggesting a finite emission time for pions. 
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I. INTRODUCTION 

Lattice calculations of the thermodynamical functions of 
quantum chromodynamics (QCD) indicate a rapid transition 
from a confined hadronic phase to a chirally symmetric de- 
confined quark gluon plasma (QGP), for a critical tempera- 
ture T c ~ 170 MeV Of. The creation and study of this QGP 
is one of the most important goals of current ultra-relativistic 
heavy ion research ||2|]. Central to this endevour, is the ques- 
tion of whether or not the QGP phase transition is first or- 
der (AT = 0), 2 nd order, or a rapid crossover reflecting an in- 
crease in the entropy density associated with the change from 
hadronic (dn) to quark and gluon (dg) degrees of freedom. 
Currently, lattice calculations seem to favor a rapid crossover 
(HIU i.e. Ar 7^ for the transition. 

An emitting system which undergoes a sharp first order 
QGP phase transition is predicted to show a larger space-time 
extent, than that for a purely hadronic system ||5|,|g]. Here, the 
rationale is that the transition to the QGP "softens" the equa- 
tion of state (ie. the sound speed c s ~ 0) in the transition re- 
gion, and this delays the expansion and considerably prolongs 
the lifetime of the system. That is, matter-flow through this 
soft region of energy densities during the expansion phase, 
will temporarily slow down and could even stall under suit- 
able conditions. 

To explore this prolonged lifetime, it has been common 
practice to measure the widths of emission source functions 
(assumed to be Gaussian) in the out- side- and long-direction 
(Rout, j^sMe and Riong) of the Bertsch-Pratt coordinate system 
iH HI S u£&. For such extractions, the Coulomb effects on 
the correlation function are usually assumed to be separable 
lHHl as well. The ratio R out /R ■side > 1, is expected for sys- 
tems which undergoe a strong first order QGP phase transition 
lHH^l]. This is clearly illustrated in Fig.Q] where the values ob- 
tained from a hydrodynamical model calculation [6] are plot- 
ted as a function of energy density (in units of T c s c ; s is the 
entropy density). These rather large ratios (especially those 
in Fig.QJ and c) have served as a major motivating factor for 
experimental searches of a prolonged lifetime at several ac- 
celerator facilities Hi 1 03 m \m H • 

An interesting aspect of Fig. UJis the strong dependence of 
Rout/Rside on the width of the transition region (i.e. AT ^ 
0). A comparison of Figs. QJ and b shows that this ratio is 



considerably reduced (by as much as a factor of four) when the 
calculations are performed for Ar = 0.1 T c . Thus, the space- 





FIG. 1: Rout I B-side as a function of the initial energy density (in units 
ofT c s c )for an expanding fireball. (a,b) are for dg/dff = 37/3, (c,d) 
for dq/df] = 3. The thick lines in (a,c) are for AT = 0, in (b,d) for 
AT = 0. 1 T c . Thin lines show results for an ideal gas case. Solid lines 
show results for T = 0.7T C , dotted for T = 0.9T C , dashed for T = T c . 
The figure is taken from Ref. J^/. 



time extent of the emitting system is rather sensitive to the 
reaction dynamics and whether or not the transition is a cross 
over. It has also been suggested that the shape of the emission 
source function can provide signals for a second order phase 
transition and whether or not particle emission occurs near to 
the critical end point in the QCD phase diagram lfl5ll . 



II. REACTION DYNAMICS AT RHIC 

Experimental evidence for the creation of locally equili- 
brated nuclear matter at unprecedented energy densities at 
RHIC, is strong [IHHIlltllli&IIlElll-Jet suppres- 
sion studies indicate that the constituents of this matter inter- 
act with unexpected strength and this matter is almost opaque 
to high energy partons 124, l25ll . 
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FIG. 2: (Color online) V2 vs. pj (left panel) and KEj (middle panel). The scaled results in the right panel is obtained via n q scaling of the data 
shown in the middle panel. Results are shown for several particle species produced in minimum bias Au+Au collisions at -vAaw = 200 GeV 
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Elliptic flow measurements 1 26, 27] validate the predictions 
of perfect fluid hydrodynamics for the scaling of the elliptic 
flow coefficient vi with eccentricity e, system size and trans- 
verse kinetic energy KEj 0281 l29l 13011 ; they also indicate the 
predictions of valence quark number {n q ) scaling |[3ll[32l[33ll . 
The result of such scaling is illustrated in Fig. |2fc; it shows 
that, when plotted as a function of the transverse kinetic en- 
ergy KEj and scaled by the number of valence quarks n q 
of a hadron (n q = 2 for mesons and n q = 3 for baryons), 
v'2 shows a universal dep endence for a broad range of par- 
ticle species 11261 1271 13411 . This has been interpreted as evi- 
dence that hydrodynamic expansion of the QGP occurs dur- 
ing a phase characterized by (i}a rather low viscosity to en- 
tropy ratio r\/s HqI |H S [33? ] and (ii) independent 
quasi-p articles which exhibit the quantum numbers of quarks 

El 111 SIM HI HI- 

The scaled V2 values shown in Fig. |2j; allow an estimate 
of c s because the magnitude of V2/E depends on the sound 
speed [30]. One such estimate [26, 2l gives c s ~ 0.35±0.05; 
a value which suggests an effective equation of state (EOS) 
which is softer than that for the high temperature QGP yfl . It 
however, does not reflect a strong first order phase transition 
in which c s ~ during an extended hadronization period. This 
sound speed is also compatible with the fact that Viipr) is ob- 
served to saturate in Au+Au collisions for the collision energy 
range y/s^ = 60 - 200 GeV 01- 

Femtoscopic measurements involving the use of the 
Bowler-Sinyukov 3D HBT analysis method [in Bertsch-Pratt 
coordinates], have been used to probe for long- rang e emis- 
sions from a possible long-lived source 

EE SUSHIS 

[3. The observed RMS-widths for each dimension of the 
emission source Ri ong ,/? s ide an d ^out, show no evidence for 
such emissions. That is, ^out/^side ~ 10. It is somewhat 
paradoxical that these observations have been interpreted as 
a femtoscopic puzzle [7] - "the HBT puzzle" - despite the 




FIG. 3: (Color online) r\/s vs (T — T c )/T c for several substances as 
indicated. The calculated values for the meson-gas have an associ- 
ated error of ~ 50% 01- The lattice QCD value T c = 170 MeV 0] 
is assumed for nuclear matter. The lines are drawn to guide the eye. 



fact that they clearly reflect a rich and complex set of thermo- 
dynamic trajectories for which c s ^ (or AT ^ 0) during an 
extended hadronization period. 

Hints for even more complicated reaction dynamics is fur- 
ther illustrated via a comparison of T[/s values for nuclear, 
atomic and molecular substances in Fig. [3] It shows the obser- 
vation that for atomic and molecular substances, the ratio r\/s 
exhibits a minimum of comparable depth for isobars passing 
in the vicinity of the liquid-gas critical point 11391 14011 . When 
an isobar passes through the critical point (as in Fig. 0, the 
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minimum forms a cusp at T c ■; when it passes below the critical 
point, the minimum is found at a temperature below T c (liquid 
side) but is accompanied by a discontinuous change across the 
phase transition. For an isobar passing above the critical point, 
a less pronounced minimum is found at a value slightly above 
T c . The value r\Js is smallest in the vicinity of T c because this 
corresponds to the most difficult condition for the transport of 
momentum i40ll . 

Given these observations, one expects a broad range of tra- 
jectories in the {T,/jb) plane for nuclear matter, to show r\/s 
minima with a possible cusp at the critical point (jjb is the 
baryon chemical potential). The exact location of this point 
is of course not known, and only coarse estimates of where it 
might lie are available. The open triangles in the figure show 
calculated values for r|/.s along the /jb = 0, rig — trajectory. 
For T < T c the T|/s values for the meson-gas show an increase 
for decreasing values of T. For T greater than T c , the lattice 
results |41] indicate an increase of r\/s with T, albeit with 
large error bars. 

These trends suggest a minimum for T|/s, in the vicinity 
of T c , whose value is close to the conjectured absolute lower 
bound ^/s = 1 /4;t. Consequently, it is tempting to speculate 
that it is the minimum expected when the hot and dense QCD 
matter produced in RHIC collisions follow decay trajectories 
which are close to the QCD critical end point (CEP). Such 
trajectories could be readily followed if the CEP acts as an at- 
tractor for thermodynamic trajectories of the decaying matter 

MM. 

The above insights on the reaction dynamics of hot and 
dense QCD matter, strongly suggest that full utility of the fem- 
toscopic probe requires detailed measurements and theoretical 
investigations of both the shape and the space-time character 
of emission source functions. Here, it should be stressed that 
a good understanding of the space-time evolution of the QGP 
is crucial. This is because, in contrast to other signals, it is de- 
rived from the influence of the equation of state on the collec- 
tive dynamical evolution of the system. This is less clear-cut 
for many other proposed signals B441, l45tl . 



III. SOURCE IMAGING: A FEMTOSCOPIC PARADIGM 
SHIFT 

A. ID source functions via imaging 

An initial strategy has been to extract the "full" distribution 
of the ID emission source function [46]. The ID Koonin-Pratt 
equation j47il : 



C(q)- l=4jt j drr 2 K (q,r)S(r), 



(1) 



relates this source function S(r) (i.e. the probability of emit- 
ting a pair of particles at a separation r in the pair center of 
mass (PCMS) frame) to the ID correlation function C(q). The 
angle-averaged kernel Ko(q,r) encodes the the final state in- 
teraction (FSI) which is given in terms of the final state wave 
function *,(r), as K Q (q,r) = ± / ^(cos(9 a( r))(|<I> q (r)| 2 - 1), 
where q r is the angle between q and r [14811 . 
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FIG. 4: (color online) (upper, a) (filled circles) Correlation function, 
C(q) for JT + 7t + and 7T~7t~ pairs; (open squares) restored correlation 
from imaging technique; (dotted line) direct correlation fitting; (open 
triangles) ID angle-averaged correlation of 3D correlation function, 
(lower) ID source function (b) S(r) and (c) 47tr 2 S(r): (open squares) 
imaging; (filled circles) spheroid fit to correlation function; (open tri- 
angles) angle-averaging of 3D-Gaussian source function. The figure 
is reproduced from Ref. 14611 



The model-independent imaging technique of Brown and 
Danielewicz j48l l49ll is used to invert Eq. [T]to obtain S(r). 
In brief, it utilizes a numerical calculation of the two particle 
wave function (including FSI) to produce an inversion matrix 
that operates on the measured correlation function to give the 
associated source function. A detailed description of the pro- 
cedure for robust measurement of the correlation functions is 
given in these proceedings ll50ll . 

The open squares in Fig.|4jb) show the ID source image ob- 
tained from the correlation function (filled circles) presented 
in Fig. EJa). The open squares in Fig. |4ja) show the cor- 
responding restored correlation function obtained via Eq. Q] 
with the extracted source function as input. Their consistency 
serves as a good cross check of the imaging procedure. 

FigureHJb)) indicates a Gaussian-like pattern at small r and 
an unresolved "tail" at large r. For comparison, the source 
function constructed from the parameters (/?iongv^side, ^out 
and X), obtained in an earlier 3D HBT analysis jl4\ is also 
shown; the associated 3D angle-averaged correlation function 
is shown in Fig. Ufa). The imaged source function exhibits a 
more prominent tail than the angle averaged 3D HBT source 
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function. These differences reflect the differences in the asso- 
ciated correlation functions for q < 10 — 15 MeV, as shown 
in Fig. Ufa). This could come from the Gaussian shape as- 
sumption employed for S(r) in the 3D HBT analysis. That is, 
the 3D Gaussian fitting procedure by construction is sensitive 
only to the main component of S(r), and thus would not be 
capable of resolving fine structure at small-q/large-r. 

The radial probabilities (4ftr 2 S(r)) are also compared in 
FigjlJc); they show that the differences are actually quite 
large. They are clearly related to the long-range contribution 
to the emission source function. 



B. ID source function via correlation function fitting 

Another approach for ID source function extraction is to 
perform direct fits to the correlation function. This involves 
the determination of a set of values for an assumed source 
function shape, which reproduce the observed correlation 
function when the resulting source function is inserted into 
Eq.Q] The filled circles in Figs.|4tb,c) show the source func- 
tion obtained from such a fit for a Spheroidal shape J48ll . 



S(r) = 



XR eS xe 4 4erfi(^) 



(8%R T R r) 



, for R >Rt, 



(2) 



where R en = l/y(l/R% — 1/1?q), Rt is the radius of the 
Spheroid in two perpendicular spatial dimensions and Rq = 
a x Rj is the radius in the third spatial dimension; a is a 
scale factor. The long axis of the Spheroid is assumed to be 
oriented in the out direction of the Bertsch-Pratt coordinate 
system. The fraction of pion pairs which contribute to the 
source X, is given by the integral of the normalized source 
function over the full range of r. The source function shown 
in Figs. |4jb) and|4jc), indicates remarkably good agreement 
with that obtained via the imaging technique, confirming a 
long-range contribution to the emission source function. 



C. Origin of the long-range contributions to the ID source 
function 



Extraction of the detailed shape of the ID source function 
is made ambiguous by the observation that several assumed 
shapes gave equally good fits which reproduced the observed 
correlation function. Nonetheless, one can ask whether or not 
a simple kinematic transformation from the longitudinal co- 
moving system (LCMS) to the PCMS can account for the ob- 
served tail? Here, the point is that instantaneous freeze-out of 
a source with R ou t/Rside ~ 1 in the LCMS would give a max- 
imum kinematic boost so that R out ~ y x R S ide in the PCMS 
(yis the Lorentz factor). Another possibility could be that the 
particle emission source is comprised of a central core and a 
halo of long-lived resonances. For such emissions, the pairing 
between particles from the core and secondary particles from 
the halo is expected to dominate the long-range emissions lf5ltl 
and possibly give rise to a tail in the source function. 



Both of these possibilities have been checked B46H : the in- 
dications are that the observations can not be explained solely 
by either one. However, further detailed 3D imaging measure- 
ments were deemed necessary to quantitatively pin down the 
origin of the tail in the source function and to determine its 
possible relationship to the reaction dynamics. 



IV. 3D SOURCE IMAGING 

In order to develop a clearer picture of the reaction dynam- 
ics across beam energies, an ambitious [ongoing] program 
was developed to measure 3D source functions in heavy ion 
collisions from AGS through CERN SPS to RHIC energies. 
This program exploits the cartesian harmonic decomposition 
technique of Danielewicz and Pratt f52l l53ll in conjunction 
with imaging and fitting, to extract the pair separation distri- 
butions in the out-, side- and long-directions. Model com- 
parisons to these distributions are expected to give invaluable 
insights on the reaction dynamics. 



Cartesian harmonic decomposition of the 3D correlation 
function 



Following Danielewicz and Pratt [52], the 3D correlation 
function can be expressed in terms of correlation moments as: 



C(q)-l=J?(q)=£ £ < 
/ a[...a ; 



(<?)<...«, (a,), (3) 



where / = 0,1,2,..., a,- = x,y,z, A l a[ a .(O q ) are cartesian 
harmonic basis elements (Q. q is solid angle in q space) and 
R l ai ai (q) are cartesian correlation moments given by 



R 



l 

(Xi ...(X/ 



(q) = 



(2/ + 1)!! [da 



4n 



^oti...a/ 



(a q )R(q). (4) 



In Eqs. [3]and|4] the coordinate axes are oriented so that z is 
parallel to the beam (i.e. the long direction), x is orthogo- 
nal to z and points along the total momentum of the pair in 
the LCMS frame (the out direction), and y points in the side 
direction orthogonal to both x and z. 

The correlation moments, for each order Z, can be calcu- 
lated from the measured 3D correlation function with Eq. [4] 
Another approach is to truncate Eq. [3] so that it includes only 
non-vanishing moments and hence, can be expressed in terms 
of independent moments. For instance, up to order I = 4, there 
are 6 independent moments: R°, R 2 x7 , R^ 2 , R^ 4 , RZ^ and ^ 7v2 - 
Here, R^ 2 , etc i s use d as a shorthand notation for R^ x , 
RL, etc. The independent moments can then be extracted as a 
function of q by fitting the truncated series to the experimental 
3D correlation function with the moments as the parameters 
of the fit. Moments have been so obtained up to order / = 8, 
depending on available statistics. The higher order moments 
(I > 4) are generally found to be small or zero especially for 
the lower collision beam energies. 

Figures[5]-H]show representative comparisons between the 
ID correlation functions, C(q), and the / = moment C(q)° = 
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R°(q) + 1 for mid-rapidity particles emitted in Au+Au and 
Pb+Pb collisions. For the AGS and SPS measurements, 
\yL — yo\ < 0.35, where yi is particle laboratory rapidity and 
yo is CM rapidity. For RHIC measurements \y\ < 0.35. The 
pr and centrality selections are as indicated in each figure. 
The 1 = moments are in good agreement with the ID cor- 
relation functions; such agreement is expected in the absence 
of significant pathology, especially those related to the angu- 
lar acceptance. Therefore, they attest to the reliability of the 
technique used for their extraction. Similarly good agreement 
was found for all other beam energies and analysis cuts. 
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FIG. 5: ID correlation function C(q) (open stars) and 1 = moment 
= R°{q) + 1 (solid stars) for low pj mid-rapidity Jl _ Jt~ pairs 
from central Au+Au collisions at yfs = 3.8 GeV. 

Correlation moments for I = 2 and / = 4 are shown for 
Au+Au collisions obtained at RHIC, in Fig. [9] The / = 2 
moments show an anti-correlation for C^ 2 (q) = R%,(q) an d 
positive correlations for C^ 2 (q) and C^(q), indicating that 
the emitting source is more extended in the out (x) direc- 
tion and less so in the side (y) and long (z) directions, when 
compared to the angle average C°(q) (also see discussion 
on imaged moments). It is noteworthy that only two of 
these I = 2 correlation moments are actually independent, i.e. 
(9) +C&( 9 ) = 0. 

The bottom panel of Fig. [9] shows that rather sizable sig- 
nals are obtained for / = 4. Here, it is C^ 2v2 correlation mo- 
ment which shows an anti-correlation. Moments were also 
obtained for / = 6. Since they were all found to be small, they 
do not have a significant influence on the shape of the source 
function. 

Figure [l0]shows correlation moments for Pb+Pb collisions 
obtained at SPS energies. The top and bottom panels show 
results for yfs = 8.9 and 17.3 GeV respectively. In each case, 
the rapidity and centrality selections are the same, but the pj 
selections are different as indicated. The open squares in the 
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FIG. 6: ID correlation function C(q) (open stars) and I = moment 
C(<?) = R (?) + 1 (solid stars) for low pj mid-rapidity pairs 
from central Pb+Pb collisions at \fs = 8.9 GeV. 



figure represent the result of a simultaneous fit to the indepen- 
dent moments with assumed shapes for the source function as 
discussed below. Within errors, the magnitudes of the higher 
order moments were found to be insignificant. 

The R\ 2 {q) an d Ryiio) correlation moments in the top pan- 
els of Fig. [I0]indicate an emission source which is more ex- 
tended in the out and long directions (compared to R°(q)). By 
contrast, the R 2 x2 and R^ (q) correlation moments in the bot- 
tom panels indicate positive correlations and hence an emis- 
sion source with big extension (compared to R° (q)) in the long 
direction. In part, this difference is related to a different kine- 
matic factor for the two pj selections. 

The correlation moments shown in Figs. [9] and [10] pro- 
vide the basis for quantitative extraction of the pair separa- 
tion distributions in the out-, side- and long-directions. This 
is achieved either by fitting the moments or by imaging them. 



B. Moment Fitting 

As discussed earlier, the 3D correlation function can be rep- 
resented as a sum of independent cartesian moments. Thus, a 
straight forward procedure for source function extraction is to 
assume a shape and perform a simultaneous fit to all of the 
measured independent moments. The open squares in the top 
panel of Fig. [10]show the results of such a fit with an assump- 
tion of a 3D Gaussian source function (termed ellipsoid since 
there is no requirement that all three radii are the same). The 
fit parameters are indicated in panel (c). Overall, the moments 
are reasonably well represented by an ellipsoid for y/s = 8.9 
GeV. However, a clear deviation is apparent for R^ 2 . A similar 
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FIG. 7: Same as Fig\6\but for yfs = 17.3 GeV. Result are averaged 
over runs with +ve and -ve magnetic fields. 
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FIG. 8: ID correlation function C(q) (open stars) and 1 = moment 
C(q)° = R°(q) + 1 (solid stars) for low pj mid-rapidity Ti~%~ pairs 
from central Au+Au collisions at y/s = 200 GeV. 



fit performed for the moments obtained for y/s =17.3 GeV 
show even larger non Gaussian deviations. 

The bottom panels of Fig. Qjjshow the results of a fit which 
assumes a linear combination of two Gaussians for the source 
function. The fit parameters are indicated in the figure. The 
two-source shape yields a much better representation of the 
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FIG. 9: 1 = 2 (top panels) and 1 = 4 (bottom panels) correlation 
moments for mid-rapidity (low pj) 7t _ 7t _ pairs from 0-30% central 
Au+Au collisions at y/s = 200 GeV. Note the change in notation; 
Ca,...a,(<?) = R ai...a,{l) 



data for y/s = 17.3 GeV. The parameters obtained from these 
and the earlier fits, can be used to generate source functions in 
the out-, side- and long-directions for comparison with mod- 
els. 



C. Moment Imaging 

A complimentary approach for obtaining the source func- 
tion distributions is that of imaging. In this case, both the 
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FIG. 10: 1 = 2 correlation moments for mid-rapidity, low pj Jt~Jt~ 
pairs from central Pb+Pb collisions for ^fs = 8.9 GeV (top panels) 
and ,/s = 17.3 GeV. The open squares represent the result of a si- 
multaneous fit of the moments with an an assumed shape for the 3D 
source function (see text). 



3D correlation function C(q) and source function S(r) are ex- 
panded in a series of cartesian harmonic basis with correlation 
moments Rr a a (q) and source moments S l a aXo) respec- 
tively. Substitution into the 3D Koonin-Pratt equation 



C(q)-1 = / drK(q,r)S(r) 



l a!...a;i 



i...a ; v 



(6) 



which relates the correlation moments R 



moments S l ai a/ 



ai.-.aM) to source 
(q). Eq. [6]is similar to the ID Koonin-Pratt 
equation (cf. Eq. [T} but pertains to moments describing dif- 
ferent ranks of angular anisotropy I. 

The mathematical structure of Eq|6] is similar to that of 
Eq. LI] Therefore, the same ID Imaging technique (de- 
scribe earlier) can be used to invert each correlation mo- 
ment R l ai a (q) to extract the corresponding source moment 



J a|...a ; 



(q). Subsequently, the total 3D source function is cal- 
culated by combining the source moments for each / as in 
equation © 
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FIG. 11: 1=0 moment C°(q) = R () (q) + 1 (top panel) and imaged 
source moment S°(r) (bottom panel) for mid-rapidity low pj 7t~Jt~ 
pairs from central Pb+Pb collisions ( y/s = 17.3 GeV). The open 
squares in the top panel represent the correlation moment restored 
from the imaged source moment HP(r) 

Figure QT| shows the imaged / = moment 5° for the in- 
put 1=0 moment C°(q) = R°(q) + 1 for mid-rapidity low pr 
7C~7C~ pairs from central Pb+Pb collisions (bottom panel). 
This source image looks essentially exponential-like. Insert- 
ing the extracted image into the ID Koonin-Pratt equation 
(Eq. [TJ yields the restored moment shown as open squares in 
the top panel. The good agreement between the input and re- 
stored correlation moments serves as a consistency check of 
the imaging procedure. 

The sources images for the / = 2 moments Ri Y and R^,. are 
v*' shown in Figs.[T2]and[T3] They show that the anti-correlation 
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FIG. 13: Same as Fig. \12\but for the 1=2 moment C xy (q) 
Open squares in top panel represent correlation moment restored 
from imaged source moment Sy y (r). 



exhibited by R\ x results in a positive I — 2 source moment con- 
tribution S xx . That is, the negative contribution of R^ x plus R° 
causes the overall correlation function to be narrower in the x 
direction compared to the angle-average correlation function. 
A narrower correlation function in x gives a broader source 
function in x; hence, the positive contribution of S^ x to S°. 

Figure Qj] shows an opposite effect for the image for the 
RL correlation moment. Here, the positive contribution of 
this correlation moment results in a negative contribution of 
the / = 2 source moment S?L The positive contribution of Ry y 
plus R° causes the overall correlation function to be broader 
in the y direction compared to the angle-average correlation 
function. The resulting narrower source function (in y) reflects 
the negative contribution of S xy to S°. 

As discussed earlier, the imaged moments are combined to 
give the resulting pair separation distributions in the out-, side- 
and long-directions. 

D. Extracted source functions 

As discussed in the previous section, 

Therefore, the net source functions in the x (out-), y (side-) 
and z (long-) directions are simply given by the sum of the 
ID source S° and the corresponding higher / contributions 5? 
where i= x, y, z. Figs. [14] - [16] show the source functions 
extracted for Pb+Pb and Au+Au collisions for the collision 
energies indicated. For the Pb+Pb collisions, source moments 
were found to be significant only for the multi-polarities / = 



and I = 2 for < p T < 70 MeV/c. For Au+Au (.2 < p T < .36 
GeV/c) significant source moments were also obtained for I — 
4. 

In Figs. [l4]and[T5]the results obtained from source imaging 
and fits to the correlation moments are shown. Clear indica- 
tions for non-Gaussian contributions to these source functions 
are apparent. A detailed analysis of these source functions 
is still ongoing. However, a number of observations can be 
drawn. First, the ratio of the RMS radii of the source functions 
in the x and y directions are 1 .3±0. 1 and and 1 .2±0. 1 respec- 
tively. This deviation from unity points to a finite pion emis- 
sion time. Here, it is important to emphasize that the rather 
low selection < pj < 70 MeV/c, ensures that the Lorentz 
factor y is essentially 1 . That is, any extension in the x direc- 
tion due to kinematic transformation from the LCMS to the 
PCMS is negligible. 

Another important observation is that the RMS pair separa- 
tion in the z direction (cf. Figs, \l4\ c) and [T5[ f)) are 1 1 fm and 
12 fm respectively. These dimensions are much larger than 
the Lorentz-contracted nuclear diameters of 3 fm at y/s = 8.9 
GeV and 1.5 fm at y/s = 17.3 GeV. In fact, they can be taken 
as the RMS pair separation resulting from the longitudinal 
spread of nuclear matter created by the passage of the two nu- 
clei. Since the latter are moving with almost the speed of light, 
one can infer a lower bound formation time for the created nu- 
clear matter. This is estimated to be ~ 8 fm/c at y/s = 8.9 GeV 
and - 10 fm/c at = 17.3 GeV. 

Figure [16] show representative source images extracted for 
Au+Au collisions at y/s — 200 GeV. For comparison, source 
functions obtained from initial calculations performed with 
the Therminator simulation code ll54ll are also plotted. The 
filled and open circles show the results from calculations with 
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FIG. 14: Tofa/ source function in (a) x (out) (b) y (side) (c) z (long) 
direction for low pj mid-rapidity 7t~7t~ pairs from central Pb+Pb 
collisions at yfs = 8.9 GeV. Source image and fits are indicated. 
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FIG. 15: Same as Fig\Mbutfor y/s = 17.3 GeV. 



and without resonance emissions turned on. Fig. [T6b shows 
that the calculations which include resonance emissions can 



account for an apparent tail in the distribution for the y (side) 
direction. However, Fig. [T6h suggest that such resonance con- 
tributions are not sufficient to account for an apparent exten- 
sion in the x (out) direction. An extensive effort to map out the 
source functions for different models and model parameters is 
currently underway. It is expected that a comparison between 
the results of such calculations and data will provide tighter 
constraints for the space-time evolution for collision systems 
spanning the AGS - SPS - RHIC bombarding energy range. 

Run4 Emc+ToF Au+Au ^s=200AGeV 




20 30 40 

r(fm) 



FIG. 16: Source images for the x (out) (a) y (side) (b) and z (long) 
( c) directions for .2 < pj < .35 GeV/c mid-rapidity Jt~ Tt~ pairs from 
central Au+Au collisions at ^/s = 200 GeV. Model calculations from 
the Therminator code are also indicated. 



V. CONCLUSIONS 

In summary, we have presented and discussed several tech- 
niques for femtoscopic measurements spanning collision en- 
ergies from the AGS to RHIC. These methods of analysis pro- 
vide a robust technique for detailed investigations of the re- 
action dynamics. Initial results from an ambitious program 
to measure a "full" femetoscopic excitation function, indicate 
that the emission source functions are decidedly non-Gaussian 
for a broad range of collision energies, and particle emission 
times can be reliably measured. 
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